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We propose an implementation of quantum neural networks using an array of quantum dots with 
dipole-dipole interactions. We demonstrate that this implementation is both feasible and versatile 
by studying it within the framework of GaAs based quantum dot qubits coupled to a reservoir of 
acoustic phonons. Using numerically exact Feynman integral calculations, we have found that the 
quantum coherence in our neural networks survive for over a hundred ps even at liquid nitrogen 
temperatures (77 K), which is three orders of magnitude higher than current implementations which 
are based on SQUID-based systems operating at temperatures in the mK range. 

PACS numbers: 03.67.Lx, 73.21.La, 72.25.Rb 


INTRODUCTION 


Quantum neural networks (QNNs) have earned 
tremendous attention recently since Google and NASA’s 
Quantum Artificial Intelligence Lab announced their use 
of D-Wave processors for machine supervised learning 
and big data classification However, these ma¬ 

chines use SQUIDs, which require temperatures in the 
mK range in order to keep quantum coherence for suffi- 
cienty long . This limitation would be lifted by orders 
of magnitude when the neural network is implemented 
with quantum dots, which we show have coherence sur¬ 
vival times on the ns scale even at 77 K. 

Quantum dots (QD) are small conductive regions of 
semiconductor heterostructure that contain a precisely 
controlled number of excess electrons, and have been the 
subject of several excellent reviews, such as The 

electrons are locked to a small region by external elec¬ 
tric and magnetic fields which define the shape and size 
of the dot, typically a few nanometers to a few hundred 
nanometers in size. Controllable QDs are often made 
on the basis of a 2D electron gas of GaAs heterostruc¬ 
tures. The energy levels of QDs are precisely controlled 
by the size of the dot and the strength of external elec¬ 
tric and magnetic fields [iEi. By arranging the quan¬ 
tum dots in a regular array on a layer of semiconductor 


heterostructure one can form a matrix for a quantum 
register, composed of either charge-based, or spin-based 
qubits, to be used for quantum computations [16|. The 
fact that QDs can easily be controlled 17| makes arrays 
of quantum dots particularly attractive for quantum neu¬ 
ral networks, where the coherence requirements are not 
as strict as in circuit-based quantum computing; instead 
the system just needs to find the minimum of energy 
state, which can be found quicker with the aid of quan¬ 
tum tunneling rather than solely classical hopping [^. 

Using an array of GaAs QDs for quantum neural net¬ 
works was first proposed by Behrman et al. [l8j |. Their 
original idea assumed the use of quantum dot molecules 
interacting with each other only by means of their shared 
phonon bath. Within this framework, it would be nearly 
impossible to control the training of the QNN since ma¬ 
nipulating a phonon bath is an arduous task [l^ . In this 
Letter we present a more achievable quantum dot based 
QNN arcitecture, where the QDs interact with each other 
via dipole-dipole coupling. We present realistic physical 
parameters for all couplings, and use a numerically exact 
Feynman integral calculation to study the time evolution 
of the phonon-damped coherence in a pair of one-electron 
QDs in such a network. 

A careful series of experiments on decoherence rates 
in SQUIDs was reported in which showed that at 
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80mK decoherence rates can be as high as 110ps“^, but 
it is well known that superconducting devices cannot op¬ 
erate at much higher temperatures. We will present nu¬ 
merically exact calculations to show that the quantum 
coherence in our quantum dot based architecture sur¬ 
vives for durations on the same time scale, but at much 
higher temperatures. 

Hypothetical QNN built of an array of QDs imple¬ 
ments the idea of a Hopfield network as the SQUID based 
processors do. To ensure potential advantages of QD 
based architecture one must prove that the three or more 
QDs interaction terms do not change dramatically the 
form of the Hamiltonian, and hence it is valid in its stan¬ 
dard Ising form H = AijSiSj + B^Si. Another required 
assumption is that the common bath of acoustic phonons 
propagating in the semiconductor heterostructure is sta¬ 
ble enough and is not significantly affected by the linear 
interaction with QDs (|1]) . Under these assumptions it is 
sufficient to show that if a pairwise interaction between 
two QDs separated by a given distance between neigh¬ 
bouring QDs allows for survival of quantum superposi¬ 
tion of states at a given temperature (of ^ lO^K order) 
we can expect the whole QD array to exhibit quantum 
coherence at high temperatures of the same order. 


INTERACTION MODEL 


We consider small QDs oi d = 3.3 nm diameter in a 
GaAs based substrate 2l|, where the excitons interact 
with their bath of acoustic phonons 21-2^. The Hamil¬ 
tonian of a pair of QD excitons embeded in a semicon¬ 
ductor heterostructure can be written in the form 


H = HE^ + Hph + V = Ho + U, (1) 

where 

2 

i=l i/j 

describes the energy of the excitons, with being the 
energy gap between the ground state and the first ex¬ 
cited state of the f-th exciton; Ki is a coupling due to an 
external driving field with Rabi frequency uil, and Jy is 
the dipole-dipole coupling between the dots, constructed 
in analogy to the dipole-dipole interaction of atoms . 
The pseudo-spin operators of the QD are 

= |X,)(X,| - |0,)(0,|, = |0,)(X,| + |X,)(0,|, 

= |X,)(0,|, = |0,)(X,|. 


The free phonon Hamiltonian is 


i?Ph = 


Pa 


m, 




^ 2mn 


( 3 ) 


and the interaction between the phonons and the QDs is 
given by 


V = Y,9axM{X,\. (4) 

cx,i 


We consider a pair of identical QDs (Ai = A 2 = 
A, J 12 = J 21 = J,Ki = K 2 = K) in which we can see 
that in the limit of vanishing driving field {K —>■ 0) the 
eigenstates of i7Ex are 


|X0) - IPX) 


|X0) -f |0X) 

V2 


|00), |XX) 


( 5 ) 


corresponding to the eigenvalues (—J, J, — A,A). The 
first two states of Eq. m have zero eigenvalue with re¬ 
spect to the interaction with phonons, given by Eq. o, 
and thus survive in coherent superposition even in the 
presence of a bath of acoustic phonons. 

After making a unitary transformation of the Hamil¬ 
tonian into a frame rotating around the z-axis with fre¬ 
quency ujl, and making a Rotating Wave Approximation, 
we arrive at a Hamiltonian in matrix form; 


/O K K 0\ 

K 0 J K 
K J 0 K ’ 

\0 K K 0/ 

where we assume K = 0.24 ps“^ is the driving field 
parameter, which corresponds to low intensity fields of 
about 0.95 kV/cm. For the dipole frequency we use the 
estimation [ 13 , [ill: J = Ki 1.4ps where 

e Ri 10 is the dielectric permittivity of GaAs, L = 7.5 nm 
is a typical inter-dot distance, and = (X0|era;|00) ~ 79 
Debye is the transition dipole moment of the QDs. 

The system of QDs interacting with phonons can be 
described in terms of the von Neumann equation for the 
reduced density matrix 

P = trph Piotij ■ ( 6 ) 

Eq. ini can be solved numerically exactly using the quasi- 
adiabatic propagator path integral technique l25ll using 
the free open source MATLAB code FeynDyn [2^. We 
use the initial condition: 


g—/3ffph 

Ptot(O) = p(0) ® (g-/3gp,) > ^here 


p{0) = |^(0))(V;(0)| , IV-(O)) = ^ (|0X) + |X0)), (7) 

and we treat the phonon bath as continuous, by defining 
the spectral density: 
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= 6{u}-uJc). ( 8 ) 

2 UlaUJa 
a 

The spectral density J{uj) for acoustic phonons in GaAs 
QDs is extremely well characterized. The form for J(w) 
has been derived from first principles, and the agreement 
with e^eriments is within the error bars of the experi¬ 
ment Q: 

J(w) = , (9) 

with a = 0.02 77 rps^ and Wc = 2.2ps“^ Q. This form 
for J(w) along with these specific parameters have been 
used consistently in a plethora of studies, in the excellent 
agreement between experiment and theory d, IH, [26l - l^ . 

In the Feynman integral representation of Eq. [BJ the 
spectral density and temperature determine the bath cor¬ 
relation function (28| : 



FIG. 1. Real part of an off-diagonal density matrix element 
with respect to time, demonstrating that the coherence sur¬ 
vives for over 100 ps even at 77 K. The graphs are shown for 
two regimes: strong coupling L = 7.5nm and moderate cou¬ 
pling L = lOnm 


R{t) = 



cos(a;t) coth 


—1 

2kBTj 


isin(a;t) , 


which in turn determines the time scale of the QD’s mem¬ 
ory. For the specific spectral density form and parameters 
used here, it has been shown that the memory lasts for 
about 2.5 ps [23,123. Therefore, having now defined all 
physical parameters used in the Hamiltonian, we deter¬ 
mined through numerical experiments that a time step of 
At = 1 ps was sufficient for numerical convergence of the 
Feynman integral, and to ensure that the full memory 
length of 2.5 ps was captured, we set the memory length 
in the Feynman integral to 3 time steps. 

Converged Feynman integral calculations for all el¬ 
ements of the QD density matrix p(t) were obtained 
at temperatures of T = 77 K and 300 K, and in Fig. 
[T] it is shown that the real part of the coherence 
Re[(X0|p(t)|X0)] lasts for over lOOps even at T = 77K. 

hile the dynamics of a coherence is expected to begin 
with a zero slope, the decoherence due to the bath causes 
the coherence to rapidly decay, in a manner closely re¬ 
sembling a decaying exponential: exp(—yt). To study 
the dependence of the decoherence rate with respect to 
temperature, we calculated the dynamics at 30 different 
temperatures between 40 K and 300 K, and for each case 
we fitted Re[(0X|p(t)|X0) ] to a decaying exponential to 
determine a decoherence rate 7 . 

It is shown in Fig. [5] that decoherence rates of 
lns“^order are still maintained even at 77K which is 
accessible by liquid nitrogen. At the same time the 
coherence dynamics essentially depends on the relative 
strength of the coupling constant (J) with respect the 
driving field (K). Changing the interdot distance, or 
varying the steepness of the QD confining potential 


Re[<0Xlp(t)IX0>] was fitted to exp(-yt) 



Temperature (K) 


FIG. 2. The decoherence rate in a pair of interacting QDs cal¬ 
culated as a linear fit of Re[{0A|p(f)|X0)] ~ exp(—yt). The 
typical decoherence rate remains at 1 ns scale. Galculated 
for d — 3.3nm, L — lOnm using 200ps overaging window 


with an external magnetic field we can achieve differ¬ 
ent regimes of coherence oscillations or suppression. The 
L = 10 nm, r = 77K curve in Fig. [1] shows the remnant 
of quantum oscillations, which is suppressed for smaller L 
(7.5 nm) or higher temperature (300K). The decoherence 
rates can be controlled in the same way. 


DISCUSSIONS 


The idea of a quantum neural network is to con¬ 
nect a set of quantum elements (in our case the QDs 
qubits) by tunable weights, so that a certain quadratic 
optimization problem can be solved. The optimization 
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then becomes a physical problem of evolving a quantnm 
system at non-zero temperature toward the minimnm en¬ 
ergy state. The dissipative bath plays an integral role in 
this quantum annealing process 30|, l3l| . This is different 
from a QNN implemented as a circuit-based quantum 
computer [^ . where the interaction with the environ¬ 
ment poses the main obstacle for creation of stable su¬ 
perpositions of quantum states. In a quantum annealer 
the interaction of the system with the environment, i.e., 
the noise, in contrast can increase the effective barrier 
transparency between the local minima and the desired 
ground state, therefore enhancing the efficiency of the 
computation 33, 3^. 

Solid state quantum annealing computers, produced by 
D-wave Systems Inc. are based on SQUID qubits 

with programmable weights implemented as inductive 
couplings. Such systems operate at the temperatures be¬ 
low IK, requiring power on the kW range for cooling the 
system. In an array of dipole-dipole coupled QDs with 
a low driving frequency the coupling weights can be 
tuned by either external fields or by changing material 
properties in the area between the dots, and we have 
shown using a numerically exact approach that such de¬ 
vices can maintain coherence at 77 K. 

The difference between our design described by the 
Hamiltonian m and the classical Hopfield neural net¬ 
work with the JijS^Sj interactions, as well as quantum 
annealers on SQUIDs, is that the interaction Ji^af a~ 
flips the states of two interacting qubits dynamically, 
in the presence of a fluctuating environment. In this 
sense our model is closer to the biological settings of the 
original Hopfield work 36| than the spin glass type en¬ 
ergy minimizing models. The Hamiltonians considered 
in this Letter can be used both for networks with self¬ 
organization and feed-forward networks [s^. 

The possible application of QDs to the construction of 
quantum neural networks has already been discussed in 
the literature [3, [l^ • These considerations however in¬ 
volve only the problem of charge control of the qubit state 
and requires manipulation of the phonon bath in order to 
work. In the present Letter, we have introduced a model 
element of a quantum dot neural network architecture. 
The QDs are assumed to be well separated from each 
other, so that wave function overlap can be neglected. 
For such array the selection rules for the interlevel tran¬ 
sitions are the same as for isolated QDs . It is possible 
to control the QDs by the near-held of plasmon modes 
propagating in the substrate beneath the dots. The QDs 
can be controlled also by external magnetic held, which 
changes the conhning potential steepness. The authors 
understand that the use of QDs with n > 2 electrons 
might be more insightful from the standpoint of using 
collective excitations of the whole QD array [iilsi , but 
such complicated models can be considered only as a fu¬ 
ture perspective for quantum neural networks [d^. We 
strongly advocate for dipole-dipole coupled QDs as an 


architecture for the construction of quantum neural net¬ 
works that are robust and feasible at higher temperatures 
than current SQUID-based architectures. 
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